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Abstract. Recent theoretical investigations (M. Beneke et al 1999 Phys. Rev. Lett. 
83 1914 ) of two-body hadronic B decays have provided justification, at least to order 
tts, for the use of the factorization ansatz to evaluate the B decay matrix elements 
provided the decay meson containing the spectator quark is not heavy. Although data 
is now available on many charmless two-body B decay channels, it is so far not very 
precise with systematic and statistical errors in total of the order of 25% or greater. 
Analyses have been made on tttt and t:K channels with somewhat contradictory results 



(M. Beneke et al 2001 Nucl. Phys. B 606 245; M. Ciuchini et al 200l|hcp-ph/01W022|). 
We take an overview of these channels and of other channels containing vector mesons. 
Because factorization involves many poorly known soft QCD parameters, and because 
of the imprecision of current data, we present simplified formulae for a wide range of 
B decays into the lowest mass pseudoscalar and vector mesons. These formulae, valid 
in the heavy quark limit, involve a reduced set of soft QCD parameters and, although 
resulting in some loss of accuracy, should still provide an adequate and transparent tool 
with which to confront data for some time to come. Finally we confront these formulae 
with data on nineteen channels. We find a plausible set of soft QCD parameters that, 
apart from three pseudoscalar vector channels, fit the branching ratios and the recently 
measured value of sin(2/3) quite well. 



PACS numbers: 12.15.Ji, 12.60.Jv, 13.25.Hw 



1. Introduction 

Much experimental effort is being expended in the study of B meson decays 0, |, 
^, and the next decade will see intensive investigation of the 5-meson system at 
the Tevatron, the SLAC and KEK i?-factories, and at the LHC The aim is to establish 
the Cabbibo-Kobayashi-Maskawa (CKM) parameters to an accuracy that will test the 
consistency of the Standard Model (SM) description of CP asymmetry, hopefully to a 
precision comparable to that of other aspects of the Standard Model. B meson decays 
should have many channels which exhibit CP asymmetry but even with precise data 
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there will be a problem of reliably unravelling the underlying weak decay mechanisms 
from the distortions caused by the strong interactions. 

Hard QCD corrections to the underlying h quark weak decay amplitudes involve 
gluon virtualities between the electroweak scale Mw and the scale 0(mb) and are 
implemented through renormalization of the calculable short distance Wilson coefficients 
Ci in the low energy effective weak Hamiltonian0 for = 1 decays at scale = 0(mb) 

HM = % I E [CMOl + C2(/i)0f] -Xt C.{fi)0, 1 + other terms (1) 

[p=u,c i=3,...,6 J 

where Xp = V*gVpb is a product of CKM matrix elements, q = d,s and the local AB = 1 
four-quark operators are 

Oi = {qa Pah-A {Pp b(3)y-A, 

03.5 = {qa ha)y-A E(^/3 ?^)vtA 

q' 

04.6 = {qa bph-A J2{^'l3 q'ah^A, (2) 

where g' G {u,d,s,c}, a and P are colour indices, and we have used the notation, for 
example, 

{qa bah-A {q'p q'lshTA = fe7/.(i - 75)M[^^7''(i T i5)q'p]- (3) 

The operators in are associated with particular processes (see for a detailed 

discussion): O12 are the tree current-current operators and Os,,,q are QCD penguin 

operators. The "other terms" indicated in (|I]) are small in the SM. They include 

magnetic dipole transition operators and electroweak penguin operators. The most 

important of these is the electroweak penguin operator 
3 

O9 = -{qa &a)v-A E ?^)v-A (4) 

where e,/ = 2/3 (—1/3) for u{d) -type quarks. The Wilson coefficient Cg is larger 
than the QCD penguin coefficients C3 and C5 and contributes —{G p / V^)XtCg{(x)OQ to 
Hes{l^)- 

Inclusion of strong interaction effects below the scale is a very difficult task 
involving, for the two-body hadronic decay B hih2, the computation of the matrix 
elements {hih2\0i\B) . Until recently, most theoretical studies of two-body hadronic 
decay invoked the factorization approximation ||^ in which final state interactions 
are neglected and {hih2\0i\B) is expressed as a product of two hadronic currents: 
(/'-il"/! /i|-B)(^2|</2'|0)- The operators O2 and 04,6 are Fierz transformed into a 
combination of colour singlet-singlet and octet-octet terms and the octet-octet terms 
then discarded. The singlet-singlet current matrix elements are then expressed in terms 
of known decay rates and form factors. Consequently, the hadronic matrix elements are 
expressed in terms of the combinations 

tt2i-l = C*2i-1 + —C2i, a2i = C*2j -|- —024-1 (5) 
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where i = 1,2,3 and = 3 is the number of colours. 



In the widely used so-called "generalized factorization" approach ||^, 0, 0, [T^, 
the renormalization scale dependence of the hadronic matrix elements {hih2\0i\B), lost 
through factorization, is compensated for through the introduction of effective Wilson 
coefficients Cf^{fi) such that 

a(/i)(o.(/i)) = cf(/i)(o.)*-^ (6) 

The effective Wilson coefficients C^^{fi),i = 3,..., 6 for the QCD penguins depend 
upon the gluon momentum and generate strong phases as crosses the uu and cc 
thresholds . The neglected octet-octet terms are compensated for by replacing Nc by 



a universal free parameter ^. The assumed universal parameters are then determined 



by fitting to as much data as possible. Some authors |jT2| have allowed the ^ parameter 
for the {V - A)iS){V - A) and {V - A) ig) {V + A) contributions to be different. 

Recently there has been significant progress in the theoretical understanding of 
hadronic decay amplitudes in the heavy quark limit [0, |15|, [16|, [1^, |1^, |19|, |20|, 0, 0, 



23| , |24| . These approaches, known as QCD (improved) factorization, exploit the fact 
that rrib is much greater than the QCD scale Aqcd and show that the hadronic matrix 
elements have the form 



{hh2\0,\B) = {h\J^^\B){h2\Jm 



l + ^r„< + 0(AQCDM) 



(7) 



provided the spectator quark does not go to a heavy meson. If the power corrections in 
Aqcd and radiative corrections in «<, are neglected, conventional or "naive" factorization 
is recovered. Although naive factorization is broken at higher order in as, these non- 
factorizable contributions can be calculated systematically. 

For B hih2 in which both hi and /12 are light, and /i2 is the meson that does not 
pick up the spectator quark, find that all non-factorizable contributions are 



real and dominated by hard gluon exchange which can be calculated perturbatively, and 
all leading order non-perturbative soft and collinear effects are confined to the B — hi 
system and can be absorbed into form factors and light cone distribution amplitudes. 
Strong rescattering phases are either perturbative or power suppressed in and, at 
leading order, arise via the Bander-Silverman-Soni mechanism [^] from the imaginary 
parts of the hard scattering kernels in ri. The contribution from the annihilation 
diagram (in which the spectator quark annihilates with one of the b decay quarks) is 
found to be power suppressed. In contrast to one of the basic assumptions of generalized 
factorization, the corrections to naive factorization are process dependent and have a 
richer structure than merely allowing ^ to be different for the {V — A) ^ (V — A) and 
{V - A) IS) {V + A) contributions. 

Similar in many ways to QCD factorization is the hard scattering or perturbative 
QCD approach |2^. However, here it is argued that all soft contributions to the 



B — hi form factor are negligible if transverse momentum with Sudakov resummation 
is included, so that the form factors can be perturbatively calculated. Also, it is argued 
that Sudakov suppression of long distance effects in the B meson is needed to control 
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higher order effects in the Beneke et al |T^, ^ approach. Non-factorizable contributions 
are now found to be complex but the imaginary parts are much smaller than that of 
the (factorizable) contribution from the annihilation diagram. Several aspects of the 
perturbative QCD approach have been strongly criticized |23|, |27[] . 

In summary, recent theoretical investigations have provided justification, at least to 
order a^, for invoking factorization to evaluate the hadronic matrix elements {hih2\0i\B) 
provided the decay meson containing the spectator quark is not heavy. Unfortunately, 
the numerical results from the QCD factorization calculations are quite sensitive to 



the input assumptions on the quark distribution functions p8[, and the contributions 
from the logarithmic end point divergences in hard spectator scatterings and weak 



annihilation p4| , thus giving rise to significant theoretical uncertainties which currently 
limit their usefulness. 



The success claimed for QCD improved factorization in fitting the observed 
branching ratios for B ^ nir and B Kit has been queried by Ciuchini et al p9[ who 
argue that, with the CKM angle 7 constrained by the unitarity triangle analysis of ||3CI|| , 
non-perturbative Aqcd/'"^6 corrections, so called "charming penguins", are important 
where the factorized amplitudes are either colour or Cabbibo suppressed. Factorization 
at leading order in Aqcd/'"^?) cannot reproduce the observed B ^ titi and B Kit 
decays and must be supplemented with enhanced c-loop penguin diagrams. 

In this paper we present formulae based on the factorization ansatz for B decay into 
the lowest mass pseudoscalar and vector mesons. We focus on the heavy meson limit 
in which combinations of the different parameters that characterize the current matrix 
elements reduce to just one. This is a consistent approach within the heavy quark 
approximation and we believe that the resulting more simple formulae will provide an 
adequate and more transparent tool with which to confront data for some time to come. 
Finally we confront these formulae with data on nineteen channels. With sin (2/3) as 
measured by the BaBar and Belle collaborations |^ we find a plausible set of soft 
QCD parameters that, apart from three pseudoscalar vector channels, fit the data well. 



2. Factorized decay amplitude 



The amplitude for charmless hadronic b decays for all these factorization schemes can 
be written as 



where 



{h,h2\H'^\B) = ^ P^mh,h2) + Q^{h,,h)] 



(8) 



Qi{hi 


h2) 


= {hi 


{Ua ba 


)v- 


-a\B) {h2\{qi3 Up)y 


-a|0) 


Q2{hi 


112) 


= {h 


{qa ba 


)v- 


-k\B) {h2\{Ui3 Ui3)y 


-a|0) 


Qsih 


h2) 


= {hi 


{qa ba 


)v- 


-a\B) {h2\{qp q'p)M~ 


a|0). 




h2) 


= {hi 


{q'a ba 


)v- 


^a\B) {h2\{qp g^)v- 


a|0), 




h2) 


= {hi 


{qa ba 


)v- 


-a\B) {h2\{qp q'p)v^ 


-a|0). 
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with 



and 



g6(/ii,M= -2(/ii|(g'„6,)s-p|B) (/i2|feg;5)s+p|0), 

Q9{hi,h2) = {hi\{qa 6„)v-a|5) (/i2|eq/(g^ '?/})v-a|0), (9) 

{q'a ba)s-P = " 75)&a, (g/3g^)s+P = ^/^(l + Ts)?^, (10) 



Pi, 2 = A„ai,2; P3,...,6,9 = •^?i«3,...,6,9 + '^cfl3,...,6,9- (H) 

The form of the matrix element Qq is a consequence of a Fierz transformation on the 
{V — A) ^{V + A) term. An example of a process in which hi (and can attach itself 
to either current is ioi q' = q = d where the matrix element for rr^p^ will have 

contributions from both terms. 

The matrix elements in Qq can be estimated from the matrix elements of the 
electroweak currents by taking the quarks, at some appropriate mass scale, to be on 
mass shell. For the current 

J^ = gi7^(lT 75)^2 (12) 

this yields 

idf.J'' = (m2 - mi)qiq2 ± (mi + m2)qil->q2- (13) 

Forming matrix elements of ([T3|) between scalar, pseudoscalar and vector states as 
appropriate, one of the terms on the right hand side will be identically zero and the 
other will then be determined by the matrix element of the left hand side. 
The Oj coefficients in (|^) have the form 

ai = + a^aP (14) 

where the leading order (LO) a}''^ are given by the naive factorization expressions 
(§) with the Wilson coefficients Cj taken at next to leading order (NLO). Detailed 



expressions for the a[^^ are given in |jT3, |T5|, |T8|, EOl E2|, E3[ for B PP, in |T9l for 



B PV and in for B VV, where P and V denote light pseudoscalar and vector 



mesons respectively. Note that claim there are errors in the expressions given by 
H, H for 5 PP. 

The beauty of the result (|14D is that the difference between the decay rate formulae 
of naive and generalized factorization as presented in numerous works 0, |12[ and 



that based upon QCD factorization lies only in the coefficients Oj, the factorization 
matrix elements Qi are common to all these approaches. Also, if the soft gluon physics 
is all accounted for in the current matrix elements then, in the heavy quark limit, Nc 
should be taken to be three. 

The corrections to the coefficients so far presented are to first order in a^. An 
encouraging feature is that they are not generally large, which leads one to hope that 
the precision with which the SM can be tested will be determined by the proximity of 
the b quark mass to the heavy quark limit and the precision of our knowledge of the soft 
QCD parameters, B meson semileptonic transition form factors and meson light cone 
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distribution amplitudes. An example of the differences in the Oj coefficients for several 
theoretical approaches is shown in Table 1. Here we compare the coefficients for QCD 
factorization, generalized factorization and a simple tree plus penguin model at scale 
Mw In this simple model, for example, the process b — > dqq produces a qq pair in a 
colour octet state so that they do not materialize in a qq meson; the q must pair with the 
d quark and the q quark with the spectator antiquark so that = = and 04 = og. 
We have estimated these coefficients by assuming a quark parton distribution function 
6x{l — x) in both hi and /i2 and that, for a given xi and X2, the qq pair has, neglecting 
quark masses, an invariant mass squared = XiX2M^. The coefficient 04 (= Og) was 
then calculated by integrating the penguin amplitude given in All estimates of 

the penguin amplitudes and the associated coefficients are uncertain; they involve 
the strong coupling ag over a distribution in g^, models of quarks in hadrons and soft 
QCD parameters, some of which are only poorly known. Because of these uncertainties, 
to be confident of any conclusions drawn about the SM it will be important to have a 
consistent picture of as many channels of charmless B decay as possible. 



3. Reduction of factorization matrix elements 



The factorization matrix elements Qi in involve products of current matrix elements 
which are evaluated through the introduction of numerous soft QCD parameters such 
as meson decay constants and transition form factors. Explicit expressions given in the 
hterature |, |Tg, for 5 ^ PP, B ^ PV and B ^ VV decay amplitudes in the 
factorization approximation are extremely cumbersome and, consequently, are unlikely 
in the short term to be of great assistance to experimentalists in analyzing their limited 
data sets. We argue that the expressions for these factorized decay amplitudes can 
be simplified, albeit at some loss of accuracy, such that a wide range of decays can 
be expressed in terms of a relatively small number of soft QCD parameters which are, 
or will be, relatively well known. Our approximations consist of (1) neglect of terms 
{mpy /rriBY in the decay amphtudes and rates and (2) use of universal generalized- 
factorization model tti coefficients rather than the (slightly) pro cess- dependent QCD- 
factorization model coefficients. The effect of this latter approximation is considered 
in section 5 for decays involving vr and K mesons for which the QCD factorization 
coefficients are well established. The loss of accuracy incurred in our approach should 
not be significant until more precise data is available. Our expressions are formally 
exact in the heavy quark limit. 

The vr and K meson decay constants are defined through the matrix elements 

(vr-|J„757'^Wa|0) = {K-\-s^^,^^'u^\Q) = -fKP^j,. (15) 

Isospin symmetry then determines that 

(vr°|M<.757'^Wa|0) = -{Adal^l^d^\Q) = -^Ul^^, (16) 



Hadronic B decays 



7 



with similar relations for the K meson. The magnitudes of /^r and fx are well determined 
from experimental measurements of the leptonic decays such as 



q2 
OTT 



mi 



(17) 



With an appropriate phase convention on the particle states, the decay constants can 
be taken to be real positive numbers which have the values /tt = (0.1307 ± 0.00046) 
GeV and fx = (0.1598 ± 0.00184) GcV. Insofar as parity is a good quantum number, 
the vector current matrix elements for tt and K are zero. By contrast, for the vector 
mesons p, u, K* and 0, the axial vector matrix elements are zero. The vector meson 
decay constants are defined by 



{u;\uc,-f''uc,\0) = {u;\da-f^dc,\0) 



{K*~\SaYUa\0) = fK*mK*e^, 

1 



V2 



(18) 



where e'^ is the meson polarization vector. Isospin symmetry determines the other matrix 
elements, for example 



(19) 



The magnitudes of some vector decay constants can be inferred from measurements of 
T lepton decay, for example 



T(t- 



K* 



' 167r •'^ ^' 



n 2 



1 + 2 



M, 



K* 



nir 



(20) 



and others from the meson decay rates into e+e pairs: 

271 .fl 



r(p° 

T{uj - 
T{(b- 



e+e") = ^a'^. 



e+e") 



e e 



rup 

a ^ 

9 3 nit^ 

2 27r 

a 

9 3 rrifj) 



(21) 



With a phase convention that makes the decay constants real and positive, these yield 
the well determined values fp = (0.216 ± 0.005) GcV, = (0.194 ± 0.004) GeV, 
= (0.233 ± 0.004) GeV and /k- = (0.216 ± 0.010) GcV. 
We now consider the B transition form factors. For pscudoscalar mesons, these are 
usually expressed in terms of two form factors Fo{t) and Fi{t): 



{p%rha\B) = 



ipB+ppY 



t 



t 



(22) 



Here = (ps — PpY ^i^d t = q^. The axial vector matrix elements are all zero. Both 
Fq and Fi are analytical functions of t with no singularity at t = 0. As there is no 
singularity in the matrix element, we have the constraint -Fo(O) = -Fi(O). The nearest 
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singularity is for t real and greater than m^, distant from t = 0. Both Fq and Fi are 
often taken as simple pole or dipole dominant. For example, from lattice QCD, 

Flit) = ^If 2^2 ^ = T^^fr-^^ (23) 

with > m^. A virtue of the rather clumsy parameterization in (^) is that when 
a contraction is taken with the matrix element for the second decay meson then only 
Fq contributes if the second meson is also a pseudoscalar and only Fi contributes if the 
second meson is a vector. For example, 

{7T-\d^^5Yu^\0){P\q/3lMB) = - UFoiml){ml-ml), (24) 
{p~Kl''u^\0){P\q(3lMB) =2\p,\mBfpF,iml) (25) 
where the momentum of the p is given by 

bp I = 



2 I 2 2 \ 2 

+ — rrip \ ^ 



2m 



B 



1/2 



(26) 



We have used e ■ pb = IPpliTiB/rrip in the B rest frame and taken the p meson to be 
moving along the z axis with zero helicity so that 

m,e^ = (pp,0,0,Ep). (27) 

If terms in {■m.„/mbY = 0.0007, {mp/mBY = 0.0212, e.t.c. are neglected, as is 
appropriate for the heavy quark limit, then, because of the analytic structure of the 
form factors and the constraint -Fo(O) = Fi{0), we can write 

{n-\ida M„)v-A|0)(P|(g/3 &/3)v-a|5) = -/,m|Fi(0) (28) 

and 

{p~\{da Uah-A\0){P\{qp &/3)v-a|5) = /pm|Fi(0). (29) 

Apart from the well determined f^^ and fp, in the heavy quark limit the B — >■ PP and 
the B PV transitions through the term ( p9D are characterized by a single parameter. 
The use of two parameters is not consistent with the heavy quark limit. 

The magnitudes of some form factors are measurable in principle through 
semileptonic decays such as 5° ^ 7r+ + + z/g, corresponding to P = 7r+, q = u, 
for which, neglecting terms proportional to the lepton mass, 

^ = if;|V'-np.l|f.WP. (30) 

However, these relations have only been used to estimate the CKM matrix element, the 
form factors have been taken from theory. For example, lattice QCD has been used to 
estimate {n~^\ua'y'^ba\B^) at large values of t where the pion is moving slowly. Various 
phenomenological forms, such as (^3|) which interpolate quite well through the calculated 
values, are then used to extrapolate to the small t region. Table 2 shows values for 
-Fvr(O) and Fk{0) from lattice QCD and other more phenomenological estimates. Other 
transition form factors are related by isospin symmetry, for example 

iAuarb^lB-) = -^(7r+|M„7'^6,|P°). (31) 
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The transition matrix elements to vector mesons are usually expressed through four 
form factors: 

(y|g,7^(l - 75)&a|5) = 2ie^,p, e'^p'^p'^s ^^^^ 



niB + my 
e^{mB + mv)Ai{t) + {pb + PvY e ■ q 
2my 



rriB + mv 



where 



+ q'" e-q 



/TUB + my 



t 



-[A,{t) - A,{t)] 



fniB - my 



A,{t). 



(32) 



(33) 



2my J \ 2my 

With an appropriate phase convention all the form factors can be taken to be real. Again 
the form factors are dimensionless analytic functions of t with the nearest singularity at 
t real and greater than m^. Also, the analytic structure demands that ^3(0) = Ao{0). 
In the semileptonic decays the matrix elements for the vector mesons to have helicity 
+1, —1 or 0, denoted by H^{t), H_{t) and -f^o(^) respectively, are given by [p^ 

2mB\py\ 



H±{t) 
Ho{t) 



{rriB + my)Ai{t) =F 



niB + my 



v{t), 



{m% - ruy - t){mB + my)Ai{t) 



2my\ft 

where the vector meson momentum \py \ in the B rest frame is 



^m\\py\^A^{t) 

VflB + Vfly 



\v\ 



' 2 I 2 J. \ 2 



2m 



B 



1/2 



(34) 



(35) 



Note that \py\i H±{t) and \^Ho(t) are analytic functions of t with singularities distant 
from t = 0. For small t it appears to be the case [B^ that Ai, A2 and V are of similar 



magnitude. Hence, for t < m^, and barring excessive cancellation, it can be anticipated 
that \/tHQ{t)/my will be larger than H±(t) by a factor of {niB/myY. 

The squares of some helicity matrix elements can in principle be measured from 
the semileptonic decays, for example, B^ —>■ + l~ + v. The decay rate for the lepton 
pair to have an invariant mass squared of t is 



[/f^(l -cos^)2 + if!(l + cos0)2 + 2/7o2sin2 9] 



(36) 



dr _ Gl t|py||\/g,p 

dt 2567r^ m\ 

where 6 is the angle between the charged lepton velocity and the recoil momentum of 
the vector boson = p+ in the lepton pair rest frame. 

As with pseudoscalar transitions, these relations (^) have only been used to 
estimate the CKM matrix element \Vub\i the form factors have been taken from theory 



such as lattice QCD. The CLEO collaboration have made such an analysis with 
several theoretical models. All models but one show a substantial dominance of i/o(^) 
at small t. Table 3 provides various theoretical estimates for the form factors and helicity 
matrix elements associated with the p+ and B^ K*^ vector decays. 
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Contraction of the transition matrix element (^) with that for a pseudoscalar 
factorization partner gives 

(P|J^|0)(F, A = 0|JJ5) = 2mB\pv\fpAo{ml). (37) 

Here A = indicates that the vector meson must have zero hehcity. For a vector 
factorization partner there are three possibihties. Since the B meson has no spin, both 
vector particles must have the same helicity so that 

{V2,\m^){ViA\J,\B) = fv,mv,H^{mlJ (38) 

where A = ±1, 0. In the heavy meson limit both mesons should have zero helicity. 
Although some cancellation between form factors can be anticipated, Table 3 suggests 
that the helicity zero states will dominate the decay rates. Also 

Mt) = (— —2 7 ) 

— rriy — t J 

so that, from (U), the constraint ^3(0) = Aq{0) and the fact that the singularities are 
distant from the small t region, we can write, for example, in the heavy meson limit 
where terms in {my /niBY are neglected, 

(vr-|(4Mjv-A|0)(r,0|(M^6/3)v-A|5) = ^m|Ao(0) (40) 

and 

{p-\{d^u^)v-A\Q){VM{u(ihp)v-A\B) = fpmlAoiO). (41) 

Thus the four soft QCD parameters characterizing decays into vector mesons reduce to 
just one in this limit. 

Based upon the simplifications discussed above, we show in tables 4 to 7 our 
expressions for the factorization matrix elements 

Qi{hi, h2) = nis^lQiihi, h2) + Qi(/i2, hi)] (42) 

for a wide range of B decays. In these tables we have used the notation Fj, = 
F[{0), Ap = Aq{0), etc. and include chiral enhancement factors = 2m^/[mfe(m„ + 
rrid)] and = 2m\ / [m[,{'mu + m^)] for the ag contributions. 



m 



B 



+ 3ml 



2my{mB + my] 



(39) 



4. Sign conventions for decay constants and form factors 

It is clear in the literature (see, for example, references |^, |10|, |32|) that different authors 
use different phase conventions for the particle states in defining the current matrix 
elements. Changes of convention should only multiply the combination I]j=i,...,6,9PiQi 
of the matrix elements Qi occurring in (Q) by a common phase factor, thus leaving 
the branching ratios unaltered. However, an inconsistent convention, such as defining 
through (vr"'"(p)|'UQ,757^(iQ,|0) = f.„p^ but insisting that /^^ is positive, will result in 
different branching ratios. If the experimental program, outlined in this paper, were 
to be carried through, the relative signs of these soft QCD parameters would not be 
determined and at this preliminary stage of B decay data analysis a theory must be 
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consulted to determine these signs. We outline here a simple quark model that illustrates 
this procedure. 

Consider the current operator j^(0) = ^7^(1 — 75)6? where the u and d quark fields 
are evaluated at = 0. To construct the matrix elements (vr"'"|j^|0) and {p'^\j'^\0) we 
take |0) to be the state with no quarks or antiquarks and |7r+) and \p~^) to be a ud pair 
at rest with a bound state S wave function (j){r) for their relative distance r. The tt"*" 
and p"*" spin states are 

k+)spin = ^,[\U, l)\d, -i) - \U, i)] (43) 

and 

|p+)spin = l)\d, -i) + \U, I)]. (44) 

We then find, up to an overall positive dimensionless factor, 

(7r+|j'^|O)cx0;(O)[-l, 0,0,0] (45) 

and 

(p+|/|O)oc0;(O)[0, 0,0,1]. (46) 
A comparison with (0), (|1^) and (p7|) then gives 

A cx 0:(O)/V^ (47) 

and 

/, (X <p;{0)/^p. (48) 

Similarly, we construct the matrix elements (7r+|'U7^(l — 75)6|5°) and (p"'"|'U7^(l — 
75)6|i?°) by assuming that the 5° is at rest with a bd S wave function ^^(r) and spin 
state 

For the 7r+ and again at rest we find, up to an overall positive factor, 

(7r+|/|5°) oc J 0;(r)0B(r)r2 dr [1, 0, 0, 0] (50) 

and 

{p-^\f\B^) cc J 0;(r)0B(r)r2 dr [0,0,0,-1] (51) 
so that, from (g^), (|32|) and (0), we infer 

Fo{{mB - m^f) oc J (f)l{r)(f)Bir)r'^ dr (52) 

and 

Ai{{mB - rrif^y) oc / (j)*{r)(f)B{r)r'^ dr. (53) 



Noting that the form factors presented in the literature do not change sign on 
extrapolation to t = then we observe that taking the wave functions to be all real and 
positive is in accord with our sign conventions for the four parameters /^r, fp, Fq and Ai. 
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5. Confronting the model with data 

The branching ratios for two-body B decays are given, in the heavy mass limit, by 

2 



(-72 3 

Bi{B^h^h2)=S,^'^^ 



327rr 



total 



Yl PiQiihi,h2] 

=1,...,6,9 



(54) 



where S = 1 unless the two bodies are identical, in which case the angular phase space is 
halved and S = 1/2. We have attempted to fit the theoretical expressions for branching 
ratios with available data. Measured branching ratios for nineteen channels are shown in 
table 8, along with references to where the data can be found. For many channels there 
are measurements presented by more than one of the three groups CLEO, BaBar and 
Belle. We have not attempted to combine the results, for each channel the table shows 
the measured branching ratio with the smallest quoted errors. We take the measured 
branching ratios to be the mean of the B and B decays: 

Br(exp) = i [Bt{B hih2) + Br(fi hih2)\ ■ (55) 

We ignore here the CP asymmetries 

^ Br(g ^ h^h^) - Br(g ^ hji^) 
CP Bi{B ^h^h2)+Bi{B ^hji2) ^ ' 

which are, so far, all consistent with being zero. 

For convenience we assign to each channel (/ii/i2) a number a. The statistical and 
systematic errors have been combined into a single error a a- The systematic errors in 
particular can be expected to be correlated. Here we ignore all correlations and form a 
function 

X^(Pj) = ^ [|Bra(Pj) — BrQ,(exp)| /cTa]^ + additional constraints. (57) 

a 

B>ia{Pi) are the theoretical branching ratios given by (|5^ ) in terms of nine parameters 
Pi,i = 1, ... ,9 which we take to be the three Wolfenstein CKM parameters {A,p,fi} 
and the six soft QCD parameters {R^ , Ptt, Fk, Ap, A^^, Ak*}- The small contribution of 
the b u penguin makes the branching ratios insensitive to i?^, which we hold fixed 
at = 0.97. The well known decay parameters {fn, fx, fp, fuj, fcf,, fx*} are held at 
their mean values and the Wolfenstein CKM A parameter is taken to be A = 0.2205. 
Additional constraints were added to the to take into account experimental and 
theoretical results lying outside the data on B decay branching ratios. For example, we 
took the Wolfenstein parameter A to be close to 0.802 as has been inferred from many 
weak decay measurements. We search for a minimum of as a function of the Pi where 
the minimum is close to the expected values of the soft QCD parameters given in tables 
2 and 3. We use the MINUIT program to minimise the x^- 

The theoretical branching ratios and contributions of the individual channels 
to x^ based on these best fit values are given in table 8. The theoretical values 
shown use the coefficients listed in table 1 as model 2. These are the process- 
independent generalized factorization coefficients computed for the renormalization 
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scale /i = mfc/2. Also in the fit we take the electroweak penguin contribution to be 
Og = —0.0094 — 0.0002i ||T^. Table 8 shows results with the systematic and statistical 
errors added in quadratureia^ = a^^^^^ + crl^sf We have also done the computation 
using simple addition of errors a = cTstat + o'syst- Both procedures are ad hoc, addition 
in quadrature reduces the influence of the systematic errors which are in general the 
smallest. The values of the best fit parameters P2,i are shown in table 9 together with 
our estimates of the two standard deviation errors. These errors are of course highly 
correlated. A plot of the error matrix ellipse for the Wolfenstein parameters p and f] 
is shown in figure 1. The results in table 9 for the best fit values of the various form 
factors lie within the spread of theoretical estimates for these form factors (see tables 2 
and 3). 



6. Discussion and conclusions 



We have investigated two-body charmless hadronic B decays within the so called QCD 
factorization model, making use of simplifications which arise from working in the heavy 
quark limit. This is particularly evident for the B —>■ VV processes which we argue to be 
predominantly of zero helicity. Consequently a wide range of decays can be expressed in 
terms of a relatively small number of soft QCD parameters, thus providing a theoretical 
framework which should be adequate to confront data for some time to come. 

Different factorization models merely modify the coefficients which premultiply 
the various combinations of soft QCD parameters, thus allowing ready comparisons 
between these models. 

If, in B decays to two vector mesons, there is a significant contribution from the 
A = ± helicity states, it should be apparent in the Dalitz type plots for the final decay 
products. Table 3 suggests that the negative helicity state might be important for B^ and 
B^ decays, and the positive helicity state for B^ and B^ decays. Since each helicity 
contributes incoherently to the branching ratio, each helicity can be considered as a 
separate channel. The additional helicity channels can be included at the cost of extra 
soft QCD parameters. The only vector channels in table 8 are B — > K*(f) and in these 
channels there is some evidence |34| for contributions from non-zero helicities. It can be 
seen from tables 1, 6 and 7 that the zero helicity amplitudes (the only ones included) 
are proportional to Ak*, a parameter which contributes significantly only to the K*(j) 
channels. Splitting the decay rates into the individual helicity channels will hardly effect 
the fit, it will only modify the estimate of Ak* given in table 9 and introduce more soft 
QCD parameters for the other helicities. 

To economize in the number of soft QCD parameters we have not included decay 
channels involving t] and f]' mesons. These amplitudes involve the mixing angle between 
the {uu + dd) and ss combinations. Also, in principle, there is mixing with cc 
which, though small, could make a significant contribution to decay modes through 
the enhanced quark decay modes h cqc. 

From table 8 it seems that of the nineteen channels included in the present analysis 
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only the PV channel 7t~K*'^, and to a lesser extent 7r~K*~^ and uK^, give a large 
contribution to the overall x^- The ujK^ channel has its largest theoretical contribution 
from the 6 — s penguin and, in particular, from the 04 and terms. The theoretical 
branching ratio is small because of the cancellation, evident in tables 1 and 6, between 
these terms. It is difficult for the theory to explain a branching ratio greater than 
2 X 10-^ 

The 7c~K*^ channel is well measured. Belle gives a large branching ratio consistent 
with the BaBar result (15.5 ± 3.4 ± 1.8) x 10"^ shown in table 8, whereas the CLEO 
value is lower. It is interesting to compare this channel with n'K^, also well measured 
and with only a marginally larger branching ratio. A fit of the theoretical ratio 

BrjB ^ TT-K^) ^ (fK \ ( ' 

^ 0.547(1 + 1.3R^Y (58) 

to the BaBar and Belle data implies < 0.5. Our fit does agree well with the CLEO 
result. 
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Table 1. Decay amplitude coefficients in the matrix elements of the effective 
Hamiltonian (||) for several theoretical models. Model 1 is a simple tree plus QCD 
penguin as described in section 2, model 2 is generalized factorization with the effective 
Wilson coefficients evaluated at the renormalization scale /i — mi,/ 2 and = ml /A 



of Beneke et al [14 



m;,/2 evaluated using the expressions 
The imaginary part of each coefficient is in parentheses. 



and model 3 is QCD factorization at /i 



iViOQcl 


„u 


„u 


u c 

Og = Og 


„u 
^4 


0.4 


1 


1.0000 


0.3333 


0.0000 


-0.0492 


-0.0524 










(-0.0234i) 


(-0.0059i) 


2 


1.0523 


0.0333 


0.0062 


-0.0397 


-0.0496 










(-0.0148i) 


(-0.0094i) 


3 


1.0510 


0.0499 


0.0050 


-0.030 


-0.038 




(0.0335i) 


(-0.1063i) 


(0.0033i) 


(-0.019i) 


(-0.009i) 




«5 = «5 










1 


0.0000 


-0.0492 


-0.0524 










(-0.0234i) 


(-0.0059i) 






2 


-0.0064 


-0.0397 


-0.0496 










(-0.0148i) 


(-0.0094i) 






3 


0.0050 


-0.050 


-0.055 








(0.0033i) 


(-0.017i) 


(-0.005i) 







Table 2. Theoretical form factors for i? ^ vr and B ^ K transitions. 



Model 



^;(0) Fk(0) 



Lattice QCD° 

Quark model wave functions'" 
Quark model wave functions'^ 
Light cone sum rule*^ 




0.27 
0.33 

0.305 



0.38 
0.341 
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Table 3. Theoretical form factors and helicity amplitudes for and 
K*^ transitions. It can be seen that, although the helicity zero channel for 
B VV decays can be expected to dominate, this table shows that much cancellation 
is anticipated in (^) . The extent to which the helicity zero states dominate should be 
apparent in Dalitz type plots for the final decay products. 



Reference 


^i(O) 


^2(0) 


V(0) 


Am 








B^- 


















1" 




0.27 


0.26 


0.35 


0.30 


0.008 


10.1 


121 


2' 




0.26 


0.22 


0.34 


0.38 


0.006 


9.4 


185 


3" 




0.27 


0.28 


0.35 


0.24 


0.008 


10.1 


81 






0.30 


0.33 


0.37 


0.21 


0.034 


11.9 


66 


5^ 




0.28 


0.28 


0.33 


0.28 


0.057 


9.9 


107 


B- - 


















2' 




0.34 


0.28 


0.46 


0.49 


0.02 


16.4 


224 






0.36 


0.40 


0.45 


0.26 


0.10 


17.1 


75 






0.37 


0.40 


0.47 


0.30 


0.30 


18.3 


93 






0.33 


0.33 


0.37 


0.32 


0.22 


13.0 


110 



Ref. ^ Ref. ^ ^ Ref. H Ref. ^ ^ Ref. ^ ^ Ref. || ^ Ref. |H| 



Table 4. Factorization matrix elements Qi{hi, ^,2) for 5° hih2 decays arising from 
b dqq. 



Decay" Qi 



Qs 



Q4 



Q9 



TT+TT 


-F f 

TT J TV 








-F f 

TV J TV 





-WF f 










F f 

TV J TV 





-F f 





-WF f 


2 


p^n~ 


-Apf-K 








Apfn 





~R^ApfTj 










2 





^pfn-\-F-jrfp 

2 





Dtt Apfn 
2 


■i{Apf^+F^fp) 
4 




-Fnfp 








F-irfp 






















Fizfuj 


DTT j4„ 




2 


2 


-^X 2 


4 


p^p- 


Apfp 








Apfp 











p^p^ 





-Apfp 





Apfp 








•^Apfp 
2 


up' 







~A.pfuj 


Aujfp+Apfuj 


~Apf^ 





fp A.pfuj 


2 


2 


4 


uouo 





A-ijfuj 






'^A^j^ 





A^f^ 
2 


K'K' 











—FxfK 





—R^FkIk 
















Ak^Ik 





—R^Ak' fx 
















FrIk* 






















AK*fK' 











07rO 








V2 





V2 





3v^ 



The decays to and 0a; are obtained from (pir' by the substitutions F^^ —>■ Ap and 
-A^ respectively. The decays to K'^K", K+K*-, K*+K- and K*+K*~ receive 
no contribution from (5i,...,6 
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Table 5. Factorization matrix elements Qdhi, /12) for B hiha decays arising from 

/; — ' (1(1(1- 

Decay" Qi Q2 Q3 Qa Qb Qe Q9 

_o_- _FVfTv_ Fyifn n n n n sFT^f-n 

^ ^ V2 V2 ^ ^ ^ ^ 2V2 

_0 „— F-^fp Apfir r\ Fi^fp — Ap^Tr r\ T31T Apf-^ 3Apf-^ 

^ P V2V2^ V2 ^ V2 2V2 

J) ^pfp ^pfp n n n n ^^pfp 

P P V2 V2 ^ ^ ^ ^ 2V2 

I ,r,- A^fp Apf^ ../o A f Auifp+Apfu, ./o A f n Apfui 

^P ^TT ~W V^^p/a; y3 V ^^pU U 

K-K'' -FKfK -R^FkIk 

K*-K^ AK*fK -R^Ak^Ik 

K-K*^ FkJk* 

K*-K*^ Ak^Ik* 

(/.TT- FJ^ F,/^ 



" The decay to 0p is obtained from that to 07r by the substitution — > Ap. 



Table 6. Factorization matrix elements Qi{h\, /12) for ^ /11/12 decays arising from 
6 — > sq'^. 



Decay Qi Q2 Q4 Qs Q 



9 



ir-7r+ -F^fK -R^FJk 

i^*-7r+ F^/i,. F,/;,. 

K-p+ ApfK ApfK 

K*-p+ ApfK* Ap/x* 

^0 <^ 

^*0^0 g ^K*f-rr g ^tt/k* g g /tt 

\/2 2\/2 

KY ^Jl^ -^0 <^ ^ 

j^*0„0 g ^K*fp g ^pIk* g g 3Aj^^^ /p 

^Oo; ^ v^F;,/^ ^ V2FkU ^ 

i^V F,,/,^ FkU FkU 

Ak*U Ak*U ^K*h 
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Table 7. Factorization matrix elements Qdhi, /12) for B hiha decays arising from 

/; sqq. 

Decay Qi Q2 Q3 Qa Q5 Qe Q9 

Tr^K*~ ^■^fK* ^K*/7r g F^rfj^^* g g SAf^* 

\/2 \/2 2'\/2 

TT-^o -R^F^fK 

7r-X*° F^/k* 

^^/i^ 

p-K*^ ApfK* 

„0l>'- ApfK Fk/p r. ^p/if n pJt: ^p/iC 3Fk/p 





Ak*U 


V2 


V2 


























>lp/K 


FkIp 


V2 


V2 




AK*fp 


V2 


V2 




FkU 


V2 


V2 




Ak*U 


V2 


V2 















K*" ^Pfn* Aii*fp g Apfif* g g 3Ajy* /p 

\/2 \/2 \/2 2\/2 

(t^K- FkI^ FkU FkU 

(t>K*- Ak'^J^ Ak^U Ak'U 
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Table 8. Measured branching ratio Br(exp), experimental error a (errors added in 
quadrature), theoretical branching ratio for best fit parameters Br(fit) and contribution 
to for various B decay channels. 



1 Idocnr 




u 




I ill I 


v2 
A 




4.1 


1.2 


Bal,Bel,Cll 


5.4 


1.07 




28.9 


6.9 


Ba2,Be2,ai 


29.6 


0.01 




3.6 


3.9 


Ba2 


0.1 


0.15 




0.8 


2.0 


Ba3,C12 


0.1 


0.15 




10.4 


3.9 


Be2,ai 


9.0 


0.13 


UJT[~ 


6.6 


2.2 


Ba3,C12 


7.15 


0.06 




10.8 


2.4 


Bal,Bel,Cll 


11.1 


0.02 




18.2 


3.9 


Bal,Bel, Cll 


18.1 


0.00 




15.9 


3.7 


Ba4,Be3,C12 


5.0 


8.75 


uoK- 


3.2 


2.2 


C12 


1.4 


0.68 


(f)K' 


7.7 


1.8 


Ba5,Be2,ai 


7.8 


0.00 


(f)K*- 


9.6 


4.4 


Ba5, Cll 


9.6 


0.00 




16.7 


2.2 


Bal, Bel, Cll 


16.9 


0.01 




22 


11 


C12 


5.0 


2.39 




8.2 


3.3 


Bal,Bel,Cll 


6.7 


0.20 




2.1 


2.1 


C12 


1.4 


0.12 




10 


6 


C12 


0.7 


2.42 




8.1 


3.2 


Ba5 


7.4 


0.05 




8.6 


3.0 


Ba5, Be2, Cll 


9.0 


0.02 



In units of 10"^. 

References to experimental groups are: 
g, Ba4 Ref. 



Ba5 Ref. §; Belle Bel Ref 



BaBar Bal Ref 

g, Be2 Ref, 



g], Ba2 Ref. [|5 



Be3 Ref. 



Ba3 Ref. 

and 



CLEO Cll Ref. 151, C12 Ref. m 



Table 9. Best fit values and one-standard deviation errors of fitting parameters. 
Fk 

0.243 ±0.031 0.312 ±0.040 0.404 ± 0.087 0.377 ± 0.067 0.349 ± 0.052 
Rl A fj p 



0.970 ±0.000 1.200 ±0.186 0.803 ± 0.063 0.375 ± 0.053 0.038 ±0.104 




Figure 1. The 95% CL error ellipse for CKM parameters p and 77 together with the 
constraints on the Unitarity Triangle from mixing, \Vub/ycb\ and sin 2/3. 



